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1^' Abstract 

c : 

Q , We show how generic off-diagonal cosmological solutions depend- 

&!)' ing, in general, on all spacetime coordinates can be constructed in mas- 

C/3 . sive gravity using the anholonomic frame deformation method. There 

^ ^ ' are found new classes of locally anisotropic and (in) homogeneous cos- 

C^ , mological metrics with open and closed spatial geometries. Such so- 

lutions describe the late time acceleration due to effective cosmolog- 
ical terms induced by nonlinear off-diagonal interactions and gravi- 
ton mass. The cosmological metrics and related Stiickelberg fields 
are constructed in explicit form up to nonholonomic frame transforms 
of the Friedmann-Lamaitre-Robertson-Walker (FLRW) coordinates. 
(^ , The solutions include matter, graviton mass and other effective sources 

QQ ' modelling nonlinear gravitational and matter fields interactions with 

^^ , polarization of physical constants and deformations of metrics, which 

may explain certain dark energy and dark matter effects. There are 
^^ , stated the conditions when such configurations mimic interesting so- 

^D ' lutions in general relativity and modifications and recast the general 

Painleve-Gullstrand and FLRW metrics. Finally, we sketch a recon- 
struction procedure for a subclass of off-diagonal cosmological solu- 
tions which describe cyclic and ekpyrotic universes. 

/\^ , Keywords: massive gravity, modified gravity, off-diagonal cosmo- 

logical solutions; ekpyrotic and little rip universe. 
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The idea that graviton may have a nontrivial mass was proposed by 
Fierz and Pauh work [Tj (for recent reviews and related f(R) modifications, 
see [2]). The key steps in elaborating a modern version of a ghost free (bi- 
metric) massive gravity theory were made in a series of papers: The so-called 
vDVZ discontinuity problem was solved using the Vainshtein mechanism [3] 
(avoiding discontinuity by going beyond the linear theory) , or following more 
recent approaches based on DGP model [4j. But none solution was found for 
another problem with ghosts because at nonlinear order in massive gravity 
appears a sixth scalar degree of freedom as a ghost (see the Boulware and 
Deser paper and similar issues related to the effective field theory approach in 
Refs. [5]). That stoppled for almost two decades the research on formulating 
a consistent theory of massive gravity. 

Recently, a substantial progress was made when de Rham and co-authors 
had shown how to eliminate the scalar mode and Hassan and Rosen estab- 
lished a complete proof for a class of bigravity / bimetric gravity theories, 
see [6]. The second metric describes an effective exotic matter related to 
massive gravitons and does not suffer from ghost instability to all orders in 
a perturbation theory and away from the decoupling limit. 

The possibility that the graviton has a nonzero mass /i results not only 
in fundamental theoretical implications but give rise to straightforward phe- 
nomenological consequences. For instance, a gravitational potential of Yuka- 
wa form ~ e~'^'~/r results in decay of gravitational interactions at scales 
r > /i~^ and this can result in the accelerated expansion of the Universe. 
This way, a theory of massive gravity provides alternatives to dark energy 
and, via effective polarizations of fundamental physical constants (in result of 
generic off-diagonal nonlinear interactions), may explain certain dark matter 
effects. Recently, various cosmological models derived for ghost free (modi- 
fied) massive gravity and bigravity theories have been elaborated and studied 
intensively (see, for instance, Refs. |3[21[S]). 

The goal of this work is to construct generic off-diagonal cosmological so- 
lutions in massive gravity theory, MGT, and state the conditions when such 
configurations are modelled equivalently in general relativity (GR). As a first 
step, we consider off-diagonal deformations of a "prime" cosmological solu- 
tion taken in general Painleve-Gullstrand (PG) form, when the Friedman- 
Lamaitre-Robertson- Walker (FLRW) can be recast for well-defined geomet- 
ric conditions. At the second step, the "target" metrics will be generated 
to possess one Killing symmetry (or other none Killing symmetries) and de- 
pend on timelike and certain (all) spacelike coordinates. In general, such 
off-diagonal solutions are with local anisotropy and inhomogeneities for ef- 
fective cosmological constants and polarizations of other physical constants 
and coefficients of cosmological metrics which can be modelled both in MGT 
and GR. Finally (the third step), we shall emphasize and speculate on im- 
portance of off-diagonal nonlinear gravitational interactions for elaborating 
cosmological scenarios when dark matter and dark energy effects can be 



explained by anisotropic polarizations of vacuum and/or de Sitter like con- 
figurations. 

We study modified massive gravity theories determined on a pseudo- 
Riemannian spacetime V with physical metric g = {g„^} and certain fiducial 
metric as we shall explain below. The action for our model is 

5 = — y 5u^^\^\[m - ^ZY(g^,,K„^) + ^L] (1) 

= -^J5u^^\^\[f{R)+ -L]. (2) 

In this formula, R is the scalar curvature for an auxiliary (canonical) con- 
nection D uniquely determined by two conditions 1) it is metric compatible, 
Dg = 0, and 2) its h- and f-torsions are zero (but there are nonzero h — v 
components of torsion 7" completely determined by g) for a conventional 
splitting N : TV =/iV©z;V, see details in [9jlj The "priority" of the con- 
nection D is that it allows to decouple the field equations in various gravity 
theories and construct exact solutions in very general forms. We shall work 
with generic off-diagonal metrics and generalized connections depending on 
all spacetime coordinates when, for instance, of type D = V + Z[7^. Such 
distortion relations from the Levi-Civita (LC) connection V are uniquely 
determined by a distorting tensor Z completely defined by T and (as a con- 
sequence for such models) by (g, N). Having constructed integral varieties 
(for instance, locally anisotropic and/or inhomogeneous cosmological ones), 
we can impose additional nonholonomic (non-integrable constraints) when 
-^IT=o ~^ ^ ^^^ R ^ R, where R is the scalar curvature of V, and it is 
possible to extract exact solutions in GR. 

The theories with actions of type ([T]) generalize the so-called modified 
f{R) gravity, see reviews and original results in [2], and the ghost-free mas- 
sive gravity (by de Rham, Gabadadze and Tolley, dRGT) [6]. We use the 
units when h = c = 1 and the Planck mass Mpi is defined via Mp^ = l/SvrG 
with 4-d Newton constant G. We write 6u'^ instead of d^u because there are 
used N-elongated differentials (see below the formulas (|lip ) and consider 
/i = const as the mass of graviton. For LC-configurations, we can fix (as a 
particular case) condition of type 

f{R)-^Uig^,,K^p) = f{R), or f{R) = f{R), or f{R) = R, (3) 

which depend on the type of models we elaborate and what classes of so- 
lutions we wont to construct. It will be possible to find solutions in ex- 
plicit form if we fix the coefficients {N^} of N and local frames for D when 



^We consider a conventional 2 + 2 splitting when coordinates are labelled in the form 
u°' = (x^,y"'), or u = {x,y), with indices i,j, k... = 1, 2 and a, b, ... — 3,4. There will be used 
boldface symbols in order to emphasize that certain geometric/physical objects and/or 
formulas are written with respect to N-adapted bases Hll[l. There will be considered also 
left up/low indices as labels for some geometric/physical objects. 



R = const and daf{R) = (dj^f) x daR = but, in general, daf{R) 7^ 0. 
The equations of motion for such modified massive gravity theory can be 
writterG 

(a^/)R^, - lf{R)g^u + ^'X^, = M-fT^,, (4) 



where Mpi is the Plank mass, R^,y is the Einstein tensor for a pseudo- 
Riemannian metric g^^ and D, T^^ is the standard matter energy-momentum 
tensor. For D — )■ V, we get R^j^ — )■ R^i, with a standard Ricci tensor i?^jy 
for V. The effective energy-momentum tensor X^j, is defined in a "sophisti- 
cate" form by the potential of graviton U = IA2 + a^U^ + a/JJi, where 03 and 
a4 are free parameters. The values U2,l^3 and U4 are certain polynomials on 



traces of some other polynomials of a matrix ICj^ = ^n ~ [ V 5 ^ ^ ) fo^' ^ 
tensor determined by four Stiickelberg fields (j)i^ as 

S^. = a^(/.^9,(/.^7?^,, (5) 

when r/^,y = (1,1,1,-1). Following a series of arguments presented in [8], 
when the parameter choice Q3 = (a — l)/3,a4 = (q^ — a + 1)/12 is useful 
for avoiding potential ghost instabilities, we can fix 

^fiu = a~^g/x!/- (6) 

De Sitter solutions for an effective cosmological constant are possible, for 
instance, for ansatz of PG type. 



ds'^ = U'^{r, t)[dr + e^/f{r, t)dtf + a'^r'^{d9^ + sin^ Od^p'^) - V'^{r, t)dt^ . (7) 
In above formula, there are used spherical coordinates labelled in the form 



u 



{x^ = r, x^ = 9,y^ = ip,y'^ = t), the function / takes non-negative 
values and the constant 5 = a/{a + 1) and e = ±1. For such bimetric 
configurations, the Stiickelberg fields are parameterized in the unitary gauge 
as (p- = t and (j)^ = rn-, (j)^ = rn-, cfA = rn-, where a three dimensional (3-d) 
unit vector is defined as n = in- = sin ^ cos 99,71- = sin0sin(^,n- = cos 6). 
Any PG metric of type ([7]) defines solutions both in GR and in MGT. It 
allows us to extract the de Sitter solution, in the absence of matter, and 
to obtain standard cosmological equations with FLRW metric, for a perfect 
fluid source 

T^,y = [p{t) + p{t)] u^Uy + p{t)g^y, (8) 

where n^ = (0, 0, 0, —V) can be reproduced for the effective cosmological 
constant ^-^-^A = ji^ /a. It is also possible to express metrics of type ([7]) in a 
familiar cosmological FLRW form (see formulas (23), (24) and (27) in [8j). 



^see details on action and variational methods in 0; we shall follow some conventions 
from [8]; the Einstein summation rule on repeating indices will be applied if the contrary 
is not stated 



Let us consider an ansatz 

ds^ = m{r,d)gi{r)dr'^ + r]2{r,0)g2{r)de^ (9) 

W(r, e, (f, t){r]3{r, 0, t)h3{r, e)[d^ + ni{r, e)dx'f 
+?74(r, 9, t)hi{r, 9, t)[dt + {wi{x'',t) + Wi{x''))dx'f}, 

with Killing symmetry on d^ = d^p, which (in general) can not be diagonalized 
by coordinate transforms. The values rja are called "polarization" functions; 
u is the so-called "vertical", v, conformal factor. The off-diagonal, N- 
coefHcients, are labelled Nf'{x^,y'^), where (for this ansatz) Nf = ni{r,6) 
and N^ = Wi{x ,t) + Wi{x ). The data for the "primary" metric are 



gi{r) = C/2 _ hi{wif,g2{r) = a\^ M = a\^sm^9,h^ = y/\fU^ - V^\, 
wi = e/7f^V^4>2 = 0,ni = 0, (10) 

when the coordinate system is such way fixed that the values /, U,V m. ([7]) 
result in a coefficient gi depending only on r. 

We shall work with respect to a class of N-adapted (dual) bases 

ea = (ej = di- N-db, Ca = da = d/dy"-) and 

e^ = {e^ = dx\e^ = dy^ + Nldy''), (11) 

which are nonholonomic (equivalently, anholonomic) because, in general, 
there are satisfied relations of type e^e^ — e^e^ = WjgGy, for certain non- 
trivial anholonomy coefficients W2a{u). For simplicity, we shall consider en- 
ergy momentum sources ([8]) and effective ([6]) which up to frame/coordinate 
transforms can be parameterized in the form 



T^ = ^(TS + a-^X;^) = ^( ™T+a-^)5g = ( ™T+ "T)5? (12) 



for constant values "^T := Mpf{djJ)'^ "T and °T = Mpf{d^J)-^a-^, 
with respect to N-adapted frames (jlip j^l 

Let us explain an important decoupling property of the gravitational 
field equations in GR and various generalizations/ modifications studied in 
details in Refs. [9j. That anholonomic frame deformation method (AFDM) 
can be applied for decoupling, and constructing solutions of the MGT field 
equations (JH with any effective source (J12p . We consider target off-diagonal 
metrics g = (^i = ??igi,^a = Vaha,]^?) [there is not summation on repeat- 
ing indices in this formula] with coefficients determined by ansatz ([9]). For 
convenience, we shall use brief denotations for partial derivatives: diip = 



^In general, such sources are not diagonal and may depend on spacetime coordinates. 
We fix such N-adapted parameterizations which will allow to construct exact solutions in 
explicit form. 



■0*,92V' = V'') f^3'0 = "0* ^-iid ^4^l^ = ip*. Computing the N-adapted coeffi- 
cients of the Ricci and Einstein tensors (see details in Refs. [9]), we transform 
([1]) into a system of nonhnear partial differential equations (PDE): 
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(t>*hl = 2/13 /i4( 
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§-^o.,., = §-,; (15) 



where =^ is used in order to show that certain equations follow from re- 
spective coefficients of the Ricci tensor R^^. In these formulas, the system 
of coordinates and polarization functions are fixed for configurations with 
91 = 92 = G and nonzero values (j)* and /i*. The equations result in solu- 

tions for the Levi-Civita configurations (with zero torsion) if the coefficients 
of metrics are subjected to the conditions 



Wi 



ej In vUmIi Sj In vl/isi = 0, diWj 



djWi and n* 



0. 



(16) 



The system of nonlinear PDE ()13p -(|16 p can be integrated in general 
forms for any to constrained by a system of linear first order equations (|14|) . 
The explicit solutions are given by quadratic elements 



ds^ 



($*)2^2 



2, ,2 



$2^^ 



( m-f + aT)$2 



4 ^ mY -j- "X) 

h4[dt + {di A)dx'f. 



hs[dip+ {dk n)dx^f 

(17) 



for any <l> = $, (9j6)* = di^* and Wi + Wi = 9jl>/6* = di Au To generate 
new solutions we can consider arbitrary nontrivial sources, ""T + "T 7^ 0, 
and generating functions, (^{x^ ,t) := e'^ and n^ = dkn{x^). Such metrics are 
generic off-diagonal and can not be diagonalized via coordinate transforms 
in a finite spacetime region because, in general, the anholonomy coefficients 
W^fl for (jlip are not zero (we can check by explicit computations). The 
polarization r/-functions for (|17p are computed in the form 



m 



-y/gi,V2 



e^/g2,m = *74 ( ''"T+ "T),ry4 = ($*)7( "T^ + "T)$^. (18) 



*We can construct exact solutions even such conditions are not satisfied, i.e. the zero 
torsion conditions are not stated or there are given in non-explicit form; this way, it is 
possible to generate off-diagonal metrics and nonholonomically induced torsions etc, see 
details in [§]. There will pe presented physical arguments for what type of generating/ 
integration functions and sources we have to chose in order to construct realistic scenarios 
for Universe acceleration and observable dark energy/ matter effects. 



So, prescribing any generating functions ^{r,6,t), n{r,6), u}{r,9,ip,t) and 
sources ™T, "T and then computing A{r,6,t), we can transform any PG 
(and FLRW) metric g = {gi,ha,Wi,hi) in MGT and/or GR into new classes 
of generic off-diagonal exact solutions depending on all spacetime coordi- 
nates. Such metrics define Einstein manifolds in GR with effective cosmo- 
logical constants '"T + "T. With respect to N-adapted frames ()lip the 
coefficients of metric encode contributions from massive gravity, determined 
by "T, and matter fields, included in ™T. 

It is possible to provide an "alternative" treatment of (J17p as exact so- 
lutions in MGT. Following such an approach, we have to define and ana- 
lyze the properties of fiducial Stiickelberg fields (j)^ and the corresponding 
bimetric structure resulting in target solutions g = (gi,ha,N'j'): Let us an- 
alyze the primary configurations related to (p^ = {(jA = a{T)pa^^n-,(/A = 
a{T) pa~^n-, (jA = tk^^), when the corresponding prime PG-metric g is 
taken in FLRW form ds^ = a^{dp^/{l - Kp^) + p'^{dd'^ + sin^ Odip'^)) - dr^. 
The related fiducial tensor ([5]) is computed 



tavdvi^u^ = ^ \dp^ + p'^(d9^ + sin^ 9dip'^) + 2HpdpdT - (^^ - H'^p'^)dT^] , 



where the coefficients and coordinates are re-defined in the form r ^ p = 
ar/a{T) and t —^ t = Kt, for K = 0, ±1; n is an integration constant; 
H := dlna/dr and the local coordinates are parameterized in the form 



X- = p,x- = 6, y- = (p,y- = r. 



For a target metric g = ga/3 and frames e„ = eirda 

for e^ 



a P 

ga/3 = e^e^rjap 



NtKb Kb 



e/ 



The values gij = e^-efqa^ = e^Sij = diag[r]igi], hab ^, ^^ .ig^p 
diag{rjaha\,Nf = d^ n and Nf = d^ A are related algebraically to data 
(|18p resulting in off-diagonal solutions ()17p . We can compute the "target" 
Stiickelberg fields as (jj^^ = e^^ (j)i^ with e'^^ being inverse to e^r, and the 

fiducial tensor Ti^p = (e^(f)!^){ei^(j)^)rifj_i, = ecTGs-'Safi- If the value S^j^ carries 
information about two constants k and a, a tensor T,^i, associated to off- 
diagonal solutions encodes data about generating and integration functions 
and via superpositions on possible Killing symmetries, on various integration 
constants [10]. In the framework of MGT, two cosmological solutions g and 
g related by nonholonomic deformationqfl are characterised respectively by 
two invariants I^ = g'^f^ da(j)-dp(j)^ and I"^ = g^^e^c/Ae^c/)^. The tensor M 
does not contain singularities because there are not coordinate singularities 
on horizon for PG metrics. The symmetry of T,^u is not the same as that 



^involving not only frame transforms but also deformation of the linear connection 
structure when at the end there are imposed additional constraints for zero torsion 



of S^^ and the singular behaviour of I"^ depends on the class of generating 
and integration functions we chose of constructing a target solution g. 

In GR and/or Einstein-Finsler gravity theories [9] , off-cosmological solu- 
tions of type pT|) were found to generalized various models of Biachi, Kasner, 
Godel and other universes. For instance, Bianchi type anisotropic cosmologi- 
cal metrics are generated if we impose corresponding Lie algebra symmetries 
on metrics. It was emphasized in [B] that "any PG-type solution in general 
relativity (with a cosmological constant) is also a solution to massive grav- 
ity. " Such a conclusion can be extended to a large class of generic off-diagonal 
cosmological solutions generated by effective cosmological constants but it 
is not true, for instance, if we consider nonholonomic deformations with 
nonholonomically induced torsion like in metric compatible Finsler theories. 

We note that the analysis of cosmological perturbations around an off- 
diagonal cosmological background is not trivial because the fiducial and refer- 
ence metrics do not respect the same symmetries. Nevertheless, fluctuations 
around de Sitter backgrounds seem to have a decoupling limit which implies 
that one can avoid potential ghost instabilities if the parameter choice is con- 
sidered both for diagonal and off-diagonal cosmological solutions, see details 
in |11] . This special choice also allows us to have a structure X^,y ~ g^^, at 
list in N-adapted frames when the massive gravity effects can be approxi- 
mated by effective cosmological constants and exact solutions in MGT which 
are also solutions in GR. 

Let us consider three examples of off-diagonal cosmological solutons with 
solitonic modifications in MGT and (with alternative interpretation) GR. 
Two and three dimensional solitonic waves are typical nonlinear wave con- 
figurations which can be used for generating spacetime metrics with Killing, 
or non-Killing, symmetries and can be characterised by additional paramet- 
ric dependencies and solitonic symmetries. 

Example 1: Taking a nonlinear radial (solitonic, with left s-label) gen- 
erating function 

$ = ^$(r, t) = 4 arctan e«'^('"-''*)+«o (19) 

and a; = 1, we construct a metric 

^^* '^^ Mr,t)[dt + {dr A)dr]'^, (20) 



I" my _|_ af^ 5,^2 



where, for simplicity, we consider n{r, 6) = 0, A{r, t) is defined as a solution 
of *$•/ ^^* = dr A and ha are given by PG-data ([10]). The generating 
function (|19p . where o"^ = (1 — u^)~^ for constants q,qQ,v, is just a 1-soliton 
solution of the sine-Gordon equation ^$** — ^^** + gin ^^ = 0. For any 



class of small polarizations with rja ~ 1), we can consider that the source 
( '"T + "T) is polarized by ^$"2 when /ig ~ h^ and h^ ~ h^i ^$*)7 ^^-^ 
with an off-diagonal term dr A resulting in a stationary solitonic universe. 
If we consider that (c?s/)~^ = ^$~^ in (|12p . we can model /-interactions 
of type ([1]) via off-diagonal interactions and "gravitational polarizations". 
In absence of matter, '"T = 0, the off-diagonal cosmology is completely 
determined by "T when ^^ transforms /i into an anisotropically polar- 
ized/variable mass of solitonic waves. Such configurations can be modelled 
if '"T^ "T. If "T^ "T, we generate cosmological models determined 
by distribution off matter fields when contributions from massive gravity are 
with small anisotropic polarization. For a class of nonholonomic constraints 
on <I> and ip (which may be not of solitonic type), when solutions ()20p are 
of type ([9]) with t/q, ~ 1 and rij, Wj ~ 0, we approximate PG-metrics of type 
([7j). Hence, by an appropriate choice of generating functions and sources, 
we can model equivalently modified gravity effects, massive gravity contri- 
butions or matter field configurations in GR and MGT interactions. For 
well defined conditions, such configurations can be studied in the framework 
of some classes of off-diagonal solutions in Einstein gravity with effective 
cosmological constants. 

Example 2: Three dimensional solitonic anistoropic waves can be gener- 
ated, for instance, if we take instead of (|19p a generating functions *<I>(r, 6, t) 
which is a solution of the Kadomtsev-Petivashvili, KdP, equations |12) . 

± «$" + («$*+ s^ *$• + g *$•••)• = 0, 

when solutions induce certain anisotropy on O]^ In the dispersionless limit 
e — >■ 0, we can consider that the solutions are independent on 9 and de- 
termined by Burgers' equation ^^* + ^^ *<!>* = 0. The solutions can be 
parameterized and treated similarly to ()20p but with, in general, a nontrivial 
term {do A)d9 after A4, when "$•/ ^$* = A* and '^^V '^* = A'. 

Example 3: Solitonic waves can be considered for a nontrivial vertical 
conformal i;-factor as in Q, for instance, of KdP type, when u = uj{r, (p,t), 
when x^ = r, x^ = 6, y^ = ip,y'^ = t, for 

± uj^^ + {dtoj + Co Cb' + eoj'")' = 0, (21) 

In the dispersionless limit e — )• 0, we can consider that the solutions are inde- 
pendent on angle ip and determined by Burgers' equation oj* +uj uj' = 0. The 
conditions (J14p impose an additional constraint eiw = uj' + 'Wi{r, 9, (p)io* + 
ni(r, 9)uj^ = 0. In the system of coordinates when Cb' = 0, we can hx W2 = 



®In a similar form, we can construct various types of vacuum gravitational 2-d and 3-d 
configurations characterized by solitonic hierarchies and related bi-Hamilton structures, 
for instance, of KdP equations with possible mixtures with solutions for 2-d and 3-d sine- 
Gordon equations etc, see details in Ref. [13| . 



and 722 = 0. For any arbitrary generating function with LC-configuration, 
$(r, 0,t), we construct exact solutions 

{dt $) 6J .^^(^^^)[^^ ^ ^Q^ A)drf, (22) 



which are generic off-diagonal and depend on all spacetime coordinates. Such 
stationary cosmological solutions are with polarizations on two angles 9 and 
93. Nevertheless, the character of anisotropics is different for metrics of type 
(|20p and (I22p . In the third class of metrics, we obtain a Killing symmetry 
on dip only in the limit cI; ^- 1, but in the first two ones, such a symmetry 
exists generically. For (|22p . the value ^ is not obligatory a solitonic one 
which can be used for additional off-diagonal modifications of solutions and 
various types of polarizations. We can provide an interpretation similar to 
that in Example 1 , if the generating and integration functions are chosen to 
satisfy the conditions 77^ ~ 1 and ni,Wi ~ 0, we approximate PG-metrics 
of type ([7]). In a particular case, we can use a conformal u-factor which 
is a 1-solitonic one, i.e. w -^ ^{r,t) = 4 arctan e'^'^'^'^"'"*^''''^'' , where o"^ = 
(1 — v'^)~^ and constants q,qo,v, defines a 1-soliton solution of the sine- 
Gordon equation w** — w" + sin a; = 0. Such a soliton propagates in time 
along the radial coordinate. 

We now consider a reconstruction mechanism with distinguished off- 
diagonal cosmological effects [9] by generalizing some methods elaborated 
for f{R) gravity in [2j. Any cosmological solution in massive, MGT and/or 
GR parameterized in a form ([9]) (in particular, as ()20p and ([22])) can be 
encoded into an effective functional / — ^U = /(i?),i?|g y = R Q. This 

allows us to work as in MGT; the conditions daf{R) = if i? = const 
simplify substantially the computations. The starting point is to consider a 
prime flat FLRW like metric 

ds^ = a^{t)[{dx'^f + {dx^f + (dy3)2] _ dt^, 

where t is the cosmological time. In order to extract a monotonically expand- 
ing and periodic cosmological scenario, we parameterize In \a{t)\ = HQt + d{t) 
for a periodic function d(t + r) = ^acos(27rt/r), where < ^a < Hq. Our 
goal is to prove that such a behavior is encoded into off-diagonal solutions 
of type dSOD-dlSD. 

We write FLRW like equations with respect to N-adapted (moving) 
frames (|lip for a generalized Hubble function H, 

SH"^ = Svrp and 3iJ^ + 2eiH = -Svrp. 

Using variables with daf{R),ji_ , = 0, we can consider a function H{t) 
when e^H = dtH = H*. The energy-density and pressure of an effective 

10 



perfect fluid are computed 



^^rHidRfrnlfiR) + ^He^dRf)) - Se^H] (23) 



i87r)-^[{dRf)-\^f{R) + 2He4{dRf) + e4e4{dRf)) + e4H] 



^-lla.^„ Jl?l , 17*1 _ /o„^-lm. l„ . /|/^^ 



In N-adapted variables, the equation of state, EoS, parameter for the effec- 
tive dark fluid is defined by 

^^_p_ / + 2H*d^f _ ^U + f{R) + 2H*d^7 
P f- 6H*d^f ^U + f{R) - 6H*d^T 

when the corresponding EoS is p = —p — {2-k)~^H* and lA{t) is computed, 
for simplicity, for a configuration of "target" Stiickelberg fields (/)^ = e^^ (j)^ 
when a found solution is finally modelled by generating functions with de- 
pendencies on t. 

Taking a generating Hubble parameter H{t) = HqI + Hi sinwt, for uj = 
2tt/t, we can recover the modified action for oscillations of off-diagonal (mas- 
sive) universe (see similar details in ^), 

/TJ 
dt[uj sm.u}t—Acos ijjt{HQ+ Hi sincjt)] exp[ifoiH sinwt]. 

(25) 
We can not invert analytically to find in explicit form R. Nevertheless, we 
can prescribe any values of constants Hq and Hi and of uj and compute 
effective dark energy and dark matter oscillating cosmology effects for any 
off-diagonal solution in massive gravity and/or effective MGT, GR. To ex- 
tract contributions of jj, we can fix, for instance, f{R) = R = R and using ([1]) 
and ([2]) we can relate f{R(t)) and respective constants to certain observable 
data in cosmology. 

The MGT theories studied in this work encode, for respective nonholo- 
nomic constraints, the ekpyrotic scenario which can be modelled similarly 
to f{R) gravity. A scalar field is introduced into usual ekpyrotic models in 
order to reproduce a cyclic universe and such a property exists if we consider 
off-diagonal solutions with massive gravity terms and/or /-modifications. 
Let us consider a prime configuration with energy-density for pressureless 
matter pm, for radiation and anisotropics we take respectively pr and p^ for 
radiation and anisotropics, n is the spatial curvature of the universe and a 
target effective energy-density p ()23p . A FLRW model can be described by 



ot72 o rPrn , Pr ^ Pa '^ , i 

?,H =8'k[-^ + — + -^- — + p\. 
a-^ a* a° a^ 
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We generate an off-diagonal/massive gravity cosmological cyclic scenario 
containing a contracting phase by solving the initial problems if u; > 1, see 
()24p . A homogeneous and isotripic spatially flat universe is obtained when 
the scale factor tends to zero and the effective /-terms (massive gravity 
and off-diagonal contributions) dominate over the rests. In such cases, the 
results are similar to those in the inflationary scenario. For recovering ([25|) . 
the ekpyrotic scenario takes place and mimic the observable universe for 
t ~ Tr/2uj in the effective EoS parameter w ^ —1 -\-smujt/3ujHi cos^ cot ^ 1. 
This allows us to conclude that in massive gravity and/or using off-diagonal 
interactions in GR cyclic universes can be reconstructed in such forms that 
the initial, flatness and/or horizon problems can be solved. 

In the diversity of off-diagonal cosmological solutions which can con- 
structed using above presented methods, there are cyclic ones with singular- 
ities of the type of big bang/ crunch behaviour. Choosing necessary types 
generating and integration functions, we can avoid singularities and elabo- 
rate models with smooth transition. Using the possibility to generate non- 
holonomically constrained /-models with equivalence to certain classes of 
solutions in massive gravity and/or off-diagonal configurations in GR, we 
can study in this context, following methods in [21 [9], , big and/or little 
rip cosmology models, when the phantom energy-density is modelled by 
off-diagonal interactions. We omit such considerations in this letter. 

In summary, we have found new cosmological off-diagonal solutions in 
massive gravity with flat, open and closed spatial geometries. We applied 
a geometric techniques for decoupling the field equations and constructing 
exact solutions in f{R) gravity, theories with nontrivial torsion and noholo- 
nomic constraints to GR and possible extensions on (co) tangent Lorentz 
bundles. A very important property of such generalized classes of solu- 
tions is that they depend, in principle, on all spacetime coordinates via 
generating and integration functions and constants. After some classes of 
solutions were constructed in general form, we can impose at the end non- 
holonomic constraints, cosmological approximations, extract configurations 
with a prescribed spacetime symmetry, consider asymptotic conditions etc. 
Thus, our solutions can be used not only for elaborating homogeneous and 
isotropic cosmological models with arbitrary spatial curvature, but also for 
study "non-spherical" collapse models of the formation of cosmic structure 
such as stars and galaxies (see also [8j). 

Deriving cosmological scenarios only for diagonalizable metrics, there are 
possibilities to discriminate the massive gravity theory form the /-gravity 
and/or GR. Mathematically, we work with certain nonlinear systems of or- 
dinary differential equations with general solutions depending on integration 
constants. Following such an approach, we positively have to modify the 
GR theory in order to explain observational data in modern cosmology and 
elaborate realistic quantum models of massive gravity. Fundamental field 
equations in gravity consists from nonlinear off-diagonal systems of PDE. 
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It is very difficult to construct off-diagonal solutions and to provide certain 
physical meaning for such spacetimes. Perhaps, such solutions are important 
at large cosmological scales because of the possibility to models polarization 
of constants and nonlinear interactions, anisotropics and generalized sources 
in modified gravity theories. A new and important feature is that the off- 
diagonal anisotropic configurations allows us to model cosmic accelerations 
and massive gravity and/or dark energy and dark matter effects as certain 
effective Einstein spaces. For a large class of such solutions, we can put the 
question: May be we do not need to consider modifications of the Einstein 
gravity but only to extend the constructions to off-diagonal solutions and 
nonholonomic systems and try to apply this in modern cosmology? This is a 
quite complicated theoretical and experimental problem and the main goal 
of this paper was to analyze such constructions from the viewpoint of mas- 
sive gravity theory when off-diagonal effects can be alternatively explained 
to other types of gravity theories. 

We also developed the reconstruction method for the massive gravity 
theory which admits and effective off-diagonal interpretation in GR and /- 
modified gravity with cyclic and ekpyrotic universe solution. The expansion 
is around the GR action if we admit a nontrivial effective torsion. For zero 
torsion constraints, it is also possible to perform off-diagonal cosmological 
models keeping the constructions in the framework of the GR theory. Our 
results indicate that theories with massive gravitons and off-diagonal inter- 
actions may lead to more complicated cyclic universe. Following such an 
approach, the ekpyrotic (little rip) scenario can be realized with no need to 
introduce an additional field (or modifying gravity) but only in terms of mas- 
sive gravity or GR. Further constructions can be related to reconstruction 
scenarios of f{R) and massive gravity theories leading to little rip universe. 
The dark energy for little rip models presents an example of non-singular 
phantom cosmology. Finally, we note that other types of non-singular super- 
accelerating universe may be also reconstructed in f{R) gravity. 
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